Doping-dependent nodal Fermi velocity in Bi-2212 revealed by high-resolution ARPES 
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The improved resolution of laser-based angle-resolved photoemission spectroscopy (ARPES) al- 
lows reliable access to fine structures in the spectrum. We present a systematic, doping-dependent 
study of a recently discovered low-energy kink in the nodal dispersion of Bi2Sr2CaCu2 08+5 (Bi- 
2212), which demonstrates the ubiquity and robustness of this kink in underdoped Bi-2212. The 
renormalization of the nodal velocity due to this kink becomes stronger with underdoping, reveal- 
ing that the nodal Fermi velocity is non-universal, in contrast to assumed phenomenology. This is 
used together with laser- ARPES measurements of the gap velocity, V2, to resolve discrepancies with 
thermal conductivity measurements. 

PACS numbers: 



In a d-wave superconductor like the high-T c cuprates, 
the electronic component of low-temperature thermo- 
dynamics is dictated by the nodes, where arbitrarily 
small excitations are permitted by the gapless nature 
of these points. An intriguing aspect of cuprate phe- 
nomenology is the so-called universal nodal Fermi ve- 
locity (vf)-OQ Along the nodal direction ((0,0)-(7r,7r)) 
the velocity measured by ARPES within 50 meV of 
Eip appears to be independent of cuprate-family or the 
number of CuC>2 layers in the compound, and is also 
nearly constant across the phase diagram - from the un- 
doped insulator, across the superconducting dome, and 
in the non-superconducting metallic state at a doping 
p>0.25 - even though other electronic properties vary 
significantly with doping. [2H1] I n addition, this univer- 
sal vp, if combined with ARPES measurement of the 
superconducting gap, leads to apparent contradiction 
with thermal conductivity observed directly in transport 
measurements, 5, 6 J suggesting that crucial information 
about the nodal quasiparticles is still missing. 

ARPES data can be represented as a convolution 
between the single-particle spectral function and the 
momentum and energy resolution of the experiment. 
Naturally, with the improved resolution of laser-based 
ARPES, the measured spectrum begins to approach the 
intrinsic spectral function, and finer structure can be re- 
vealed. Recent laser-ARPES measurements along the 
nodal direction of optimally-doped Bi-2212 have uncov- 
ered a low-energy (<10 meV) kink,|7] in addition to 
the larger kink seen at 50-80 meV in all cuprates. [MTU] 
Other laser-ARPES studies have shown a correspond- 
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FIG. 1: (a)-(c) False-color image plots of the nodal dispersion 
of UD55, UD65, and UD92, measured at 10K. Solid curves 
indicate band dispersions derived from MDC peak positions. 
Inset of (a): Brillouin zone schematic with Fermi surface in- 
dicated in red. The a-axes in (a)-(c) correspond to momen- 
tum along diagonal blue line ('nodal cut') in inset, (d) EDCs 
(solid) and symmetrized EDCs (dashed) at k^. A single peak 
at Ef in the latter confirms that spectra in (a)-(c) are un- 
gapped. kF determined from Fermi crossing of dispersion 
(vertical dashed lines in (a)-(c)). 



ing decrease in nodal linewidth at low energies. [TTl - fT5] 
In this letter, we present the systematics of the low- 
energy kink by means of a doping-dependent study of 
underdoped Bi-2212. 20 samples with 6 dopings in the 
range 0.076<p<0.14 were measured using a 7 eV laser 
and a Scienta SES2002 analyzer. 7 eV photons were 
produced by second harmonic generation from a 355 nm 
laser (Paladin, Coherent, Inc.) using a nonlinear crystal 
KBe2B03F2.[14] Energy and momentum resolution were 
3 meV and better than 0.005 A -1 , respectively. Sam- 
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FIG. 2: (a) Band dispersions from Fig [T] offset for clarity. 
Dashed line accompanying UD55 denotes assumed linear bare 
band. Black dotted lines are linear fits 30-40 meV (v m id), 
extrapolated to Ef. This differs from vp (fit 0-7 meV), in- 
dicated on UD65 dispersion by orange dashed line. Pink bar 
marks 70 meV kink, and grey bar marks low-energy kink, 
where vp deviates from v m id- Inset: UD92 dispersion below 
20 meV. (b) ReE, approximated by subtracting a linear bare 
band from dispersions in (a) , is peaked at the position of the 
70 meV kink and changes slope near the energy position of 
the low-energy kink, (c) Detail of the low-energy portion of 
ReE. Thick dotted lines are fits near Ef- All curves devi- 
ate from these lines between 6-10 meV, as highlighted by the 
shaded bar. For the two most underdoped samples, the slope 
of ReE evolves between the low-energy kink and 20-30 meV, 
suggesting an additional kink, (d) 2ImE, with dotted lines as 
guides-to-the-eye. All curves decrease more rapidly near Ef, 
but most markedly in UD92. For more underdoped samples, 
the effect may be obscured by a larger linewidth and possible 
additional kink. 



pies were cleaved in-situ at a pressure <4xl0 11 torr to 
obtain a clean surface, and measured at 10K. 

Fig. [IJa)-(c) show ARPES image plots for nodal cuts 
at three dopings: UD55 (underdoped, T C =55K), UD65, 
and UD92, corresponding to hole-dopings of approxi- 
mately 0.088, 0.10, and 0.14. Standard momentum dis- 
tribution curve (MDC) analysis-Lorentzian fits at fixed 
energy- is used to extract the band dispersions. [TS] The 
energy distribution curves (EDCs)-intensity as a func- 
tion of energy at fixed momentum- at kf in Fig. [ljd) 
indicate that these spectra are ungapped, as the sym- 
metrized EDCs [16] have a single peak at Ef . 

The systematics of the low-energy kink are studied via 
the MDC-derived nodal dispersion, which are plotted for 
three dopings in Fig. [2]ja). In addition to the large, 



ubiquitous kink near 70 meV, a smaller kink is also evi- 
dent: the dispersion within 10 meV of Ef deviates from 
the velocity fit between 30-40 meV. This deviation ap- 
pears more pronounced for more underdoped samples. 
Consistent with the work of Plumb et al, the velocity 
(slope of the MDC dispersion) within 7 meV of Ef is 
smaller than the velocity at higher binding energy, no- 
tably opposite to the expected effects of instrument and 
thermal broadening. [7J We also note that the low-energy 
kink cannot be identified as an artifact due to a gap, be- 
cause measurements are performed at the node where the 
superconducting gap is zero. 

Another way to visualize the low-energy kink is via 
the real part of the electronic self-energy, ReS, plot- 
ted in Fig. U»-(c). The low -energy kink is marked 
by a deviation of the slope of ReE at 6-10 meV from 
the slope established at Ej?. For the sample closest to 
optimal doping, UD92, there is a single 'knee' in ReE. 
Meanwhile, the slope of ReE for UD55 and UD65 con- 
tinues to evolve until 20-30 meV, possibly suggesting an 
additional kink, reminiscent of the 70meV kink, which 
may have several components. |17H19| From the Kramers- 
Kronig relation between ReE and ImE, a signature of the 
low-energy kink is expected in ImE, which is proportional 
to the MDC FWHM. In Fig. gd) we show that all dop- 
ings exhibit a downturn in ImE near Ef , though this is 
most pronounced for UD92. For more underdoped sam- 
ples, the larger linewidth and possible additional kink 
make it more difficult to get quantitative information 
from ImE, but the observation that ImE decreases more 
rapidly close to Ef remains robust. The appearance of 
a low-energy feature in both ReE and ImE strongly ar- 
gues against a spurious origin for the low-energy kink, 
and the phenomenology reported in Fig. [2] is reproduced 
in the other samples in our study. Thus, the systematics 
of a new energy scale can be added to the hierarchy of 
multiple energy scales in the cuprates.|20j 

The ubiquity of the low-energy kink in UD Bi-2212 
leads us to reexamine previous measurements of vp, as 
there is now compelling evidence that quasiparticles very 
close to Ef experience a heretofore unconsidered mass 
renormalization. The nodal Vf is plotted in Fig. [3j and 
our key finding is that vp is not universal, but rather, 
has a pronounced doping dependence in the regime of 
this study. To characterize our data, at least two ve- 



locities are needed: 



the linear fit between 30-40 



meV, and vp , the velocity fit between 0-7 meV, as de- 
fined in Fig. [5Ja) . These energy ranges are chosen to get 
sufficient data points while avoiding the low-energy kink 
and 70 meV kink. v m id is found to be approximately 
1.8 eVA, without a distinct doping dependence, con- 
sistent with the previously reported 'universal' value. [1] 
Meanwhile, vp decreases monotonically with underdop- 
ing. This is consistent with recent quantum oscillation 
results, which suggest a divergence of the effective mass 
in the underdoped regime. [2"T] The coupling strength of 
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FIG. 3: (a) Doping dependence of vp and Vmid. Open 
red (blue) squares denote average Vjv (v m id) f° r each dop- 
ing. Boundary of shaded region denotes doping dependence of 
T c . Doping is approximated from empirical relation T c =96(l- 



82.6(p-0.16) i 



id has little systematic doping depen- 



dence, while vf decreases with underdoping. (b) Doping de- 
pendence of v m id/vF, which is related to the renormalization 
coupling strength. 



this low-energy renormalization can be roughly assessed 
by the velocity ratio v m id/vF, which is plotted in Fig. 
[3]jb), and suggests that coupling strength increases with 
underdoping. Notably, the ratio of velocities on either 
side of the 70meV kink exhibits the same doping depen- 
dence, though with the 70 meV kink, it is the higher en- 
ergy velocity (cj>70meV) which is doping-dependent. [1 
Although a doping-dependent vf presents a significant 
shift from our previous understanding of cuprate nodal 
physics, our results are not inconsistent with previous 
measurements: v m id is indeed doping-independent, and 
inferior energy resolution can easily obscure subtle low- 
energy kinks near Ep. Our finding underscores the im- 
portance of very low energy scales in these systems and 
revises cuprate phenomenology by linking nodal vp to 
doping and T c , previously suggested by the temperature 
dependence of the low-energy kink. [7] Further, this dop- 
ing dependence constrains the origin of the low-energy 
kink, and may aid interpretation of bulk thermodynamic 
measurements, particularly thermal conductivity, which 
will be the focus of the remaining discussion. 

For the cuprates, thermal conductivity near T=0 can 
be expressed in terms of two components of the Fermi 
velocity: the velocity perpendicular to (vp) and tan- 
gential to (v 2 ) the Fermi surface (FS) at the node (Fig. 
[4^ a)-fc)). [2"3"] - f2"5] For a 2D rf-wave superconductor in the 
clean limit, the residual linear term (T=0 extrapola- 
tion) of thermal conductivity, Kq/T, is independent of 
the quasiparticle scattering rate, interaction energy, or 
other sample-dependent parameters. [24] [26] 

In this regime, k /T is related to vp and v 2 by a simple 
formula: [Ml 121 IMI 



Kq/T 



V 2 



k\ n .Vp 
3h d v 2 vp 



(1) 



where n is the number of Cu0 2 planes per unit cell, 
and d is the c-axis unit cell length. The second term 
is usually negligible, as v 2 <§;vi?. Thus, by measuring 
bulk thermal conductivity, one can extract a microscopic 
parameter, v 2 /vp, which fully determines the ground 
state nodal electronic structure of cuprates. ^5] Ther- 
mal conductivity measurements on Lar>-^Sr^CuOA. |27l 
EH| YBa 2 Cu 3 7 ,[28] Bi 2 Si 2 - x ha x C\iO G+ s,^ and Bi- 
2212[5] have all shown that Ko/T decreases with un- 
derdoping, implying that vp/v 2 also decreases. How- 
ever, using v 2 reported in Ref. I30j and a universal vp, 
Sun et al argued that ARPES suggested a different dop- 
ing dependence of vf/v 2 . [5] Such contradictions have 
been attributed to disorder effects, such as electronic in- 
homogeneity [5] or disorder-induced magnetism, [5] but 
previous analysis lacked a crucial component: a doping- 
dependent Vp. 

For comparisons via Eqn. (1), we have obtained v 2 
from laser- ARPES measurements of the momentum de- 
pendence of the superconducting gap near the node (Fig. 
|4]Jd)), and these values are consistent with recently pub- 
lished data. [3"Tl I5"2"] Using these v 2 together with vp 
from Fig. [3]ja), we plot the ratio vp/v 2 in Fig[4|f) along- 
side the thermal conductivity values reported by Sun et 
aZ.[5] The ARPES vp/v 2 decrease strongly with under- 
doping, exhibiting a consistent trend with the thermal 
conductivity results for Bi-2212 as well as other cuprates. 
The vf/v 2 derived from ARPES and thermal conductiv- 
ity differ in absolute value, and this may be related to 
the in-plane anisotropy of thermal conductivity that has 
been reported in Bi-2212 (n a ^Kb). [33] The doping- 
dependent data in Ref. [5] were along the a-axis, which 
has a larger Ko/T, but the discrepancy in absolute value 
of vp / v 2 shown in Fig. [4] (f ) may indicate that the 6-axis 
may be the correct one to compare to ARPES. Notably 
ARPES does not observe in-plane anisotropy of nodal 
single particle parameters (vp and v 2 ), so the anisotropy 
in thermal conductivity must have a different origin, such 
as coexisting density- wave order. [34j Alternately, the dif- 
ference in absolute value olvp/v 2 may suggest a propor- 
tionality constant of different origin is needed in Eqn. (1) . 
Nevertheless, with the doping-dependent v F elucidated 
by laser ARPES, we are able to reproduce the doping- 
dependence of Ko/T which is seen in multiple cuprate 
families. 

The superior resolution of laser ARPES allows access 
to some of the lowest energy scales of high-T c cuprates, 
and the findings may further constrain the microscopic 
theory of high-T c superconductivity. With this pow- 
erful probe, we have been able to uncover the doping 
dependence of a microscopic parameter, vp. A non- 
universal nodal Fermi velocity both answers old questions 
and introduces new ones: discrepancies between doping- 
dependent ARPES and thermal conductivity measure- 
ments can be resolved, but there are new questions about 
the origin of the low-energy kink, the explanation of its 
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FIG. 4: (a) Cutout showing FS (top) and dispersion perpendicular (v>, left) and tangential (v2, right) to FS at node, from 
measurement on UD92 at 10K. (b)-(c) Image plots showing measured vj? and V2 directly. The latter image consists of EDCs 
at Icf from many parallel cuts near the node, (d) Comparison of synchrotron- and laser-based ARPES measurements of the 
superconducting gap of UD92 around the FS. When the data in (d) is fit to a simple d-wave form, A(6*)=Aocos(2#) close to 
the node, V2 ~ 2Ao/kF, where ki? is the distance from the node to (7r,7r). (e) V2 from synchrotron- and laser-based ARPES 
experiments, (f) Comparison between vf/v2 from laser ARPES and thermal conductivity from Ref. [5 . Doping dependence 
is consistent, though absolute values differ. Open squares indicate average laser ARPES values. 



doping and temperature dependence, and its role in su- 
perconductivity. 
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